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Abstract 

We define a class of pseudo-ergodic non-self-adjoint Schrodinger opera- 
tors acting in spaces and prove some general theorems about their 
spectral properties. We then apply these to study the spectrum of a non- 
self-adjoint Anderson model acting on /'^(Z), and find the precise condition 
for to lie in the spectrum of the operator. We also introduce the notion of 
localized spectrum for such operators. 
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1 Introduction 



Recent papers have obtained some striking results concerning the spectral prop- 
erties of the non-self-adjoint (nsa) Anderson model, which models the growth of 
bacteria in an inhomogeneous environment, |10, 11, 12, |l^, M. To be more precise 



the authors have determined the asymptotic limit of the spectrum of a nsa random 
finite periodic chain almost surely as the length of the chain increases to infinity. 
In a later paper the author considered the same random operator H acting on 
/^(Z), and found that the spectrum is very different from that obtained by the 
cited authors, |P. The reason for this is that the spectral properties of nsa op- 
erators are highly unstable, and infinite volume limits should be examined using 
pseudospectral ideas, [|l|, ||, ||, ||, |^, |l^ |15|, |16|, 0]. More specifically if A lies 
in the spectrum of the infinite volume nsa Anderson model, it need not be close 
to the spectrum of the finite volume periodic Anderson model; one expects rather 
that the norm of the resolvent operator [H — A)^^ of the finite volume model will 



diverge as the volume increases. These pseudospectral ideas have been worked out 
in detail for a random bidiagonal model, which is in a certain sense exactly soluble, 



1^, ^, |T8[. Our results may therefore be interpreted as finding the region in the 
complex plane for which the finite volume nsa periodic Anderson model has very 
large resolvent norm. 

In the present paper we reconsider such problems in a more general context, in 
which the probabilistic aspects have been eliminated in favour of what we call 
pseudo-ergodic ideas. As well as making the subject more accessible to those 
without a probabilistic training, this emphasizes the fact that the spectral matters 
which we consider depend only on the support of the relevant probability measure. 
On the other hand the asymptotics of the spectrum of the finite volume periodic 
nsa Anderson model does depend on the probability measure. We finally carry 
out a more detailed spectral analysis of the infinite volume nsa Anderson operator, 
and find precise conditions under which zero almost surely lies in the spectrum. 
We also obtain further results on the location of the spectrum, which come close 
to a complete determination in many cases. In the final section we consider the 
possibility that there may be constraints on the pair of values of the potential at 
two neighbouring points which are absolute rather than just probabilistic. 



2 The general context 

The operators which we consider act on the Hilbert space 1^{X,)C) ~ ® /C, 

where X is a countable set on which a group F acts by permutations. The simplest 
choice of the auxiliary Hilbert space /C is C, but other choices are needed in some 
applications; see the end of Section 3. Many of the results presented here apply 
to 1^{X, /C) with p ^ 2 without modification (the case p = 1 is of probabilistic 
importance), but this does not apply to those involving numerical ranges. We 
define the unitary operators for 7 G F by U^f{x) = /(7~^x). The bounded 
operators which we study are of the form H = Hq ^ I + V. Here Hq acts on l'^{X) 
and commutes with the action of F in the sense that HqU^ = U^Hq for all 7 G F, 
or equivalently 

Ho{ix,-iy) = Ho{x,y) 

for all 7 G F and all x,y & X, where Hq{x, y) is the infinite matrix associated with 
Hq. We assume that the spectrum E of Hq is known. From this point onwards we 
write -f^o ^01 Hq ® I. 

Given a norm closed, bounded set M. C C{}C), we assume that the operator V is 
of the form 

iVf)ix) = Vix)fix) 

where V{x) G Ai for all x E X. We say that V is (F, A^) pseudo-ergodic if its set 
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of spatial translates is dense in the following sense. For every £ > 0, every finite 
subset F C X and every W : F ^ A4, there exists 7 e F such that 

\\W{x) -VMW <e 

for all X E F. It is well known that a large class of suitably defined random 
potentials have this property almost surely, but we consider a single potential, and 
do not need to introduce any probabilistic ideas. The same class of pseudo-ergodic 
potentials is applicable to a variety of different random models, as we explain in 
more detail in the final section. 

The above definition suffices for our purposes, but it does not capture the full sense 
of random behaviour and may be refined as follows. We define a direction U to be 
an infinite subset of X such that for every finite F d X there exists 7 G F such 
that 7F C U . We then say that V is (F, M) pseudo-ergodic in the direction U if 
for every £ > 0, every finite subset F <Z X and every W : F ^ M., there exists 
7 e F such that 7F C U and 

\\W{x)-V{^x)\\ <e 

for all X & F. Suitably defined random potentials have this property for every 
choice of direction almost surely, and therefore have the property simultaneously 
for any countable set of directions almost surely. The property itself, however, is 
defined for a single potential and makes no mention of probability. 

The following theorem is an adaptation of a well-known result of Pastur for random 
potentials. We will use it to approximate Spec{H) from inside by making suitable 
choices of W. 

Theorem 1 If H = Hq + V where V is (F, A4) pseudo-ergodic and K = Hq + W 
where W : X ^ A4 is arbitrary, then 

Spec(X) C Spec(H). 

In particular ifV, W are both (F, A4) pseudo-ergodic then they have the same spec- 
trum. 

Proof If A G Spec(i^') then there exists a sequence /„ G P{X,IC) with ||/„|| = 1 
and either \\Kfn — \fn\\ — >■ or \\K*fn — Xfn\\ — >■ 0; we consider only the former 
case, the latter being similar. Given £ > a truncation procedure shows that there 
exists / with finite support F in X such that ||/|| = 1 and — A/|| < e/2. Since 
V is pseudo-ergodic there exists 7 G F such that \\Hjf — Kf\\ < e/2, where 

H^ = U-^HU^ = Ho + V{-f). 
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Putting — Ujf we deduce that 

||if/,-A/,|| = \\U;'HU,f-\f\\<e 
and the arbitrariness of £ > imphes that A e Spec{H). 
CorollEiry 2 If H — Hq + V where V is (T^M.) pseudo-ergodic then 

Spec(iJ) = |J{Spec(ifo + W):W e M^}. 
If also H — Hq + V where V is (F, Ai) pseudo-ergodic with M. C. Jv[ then 

Spec(//) C Spec(^). 

From this point we assume that H = Hq -\-V where V is (F, A^) pseudo-ergodic. 
We put 

Spec(A<) = U Spec(A) 

and 

Num(A^) = y Num(A) 

where Num denotes the closure of the numerical range. 

Theorem 3 The spectrum of H satisfies 

E + Spec(A^) C Spec(i/) C Num(i/o) + Conv(Num(Al)) 

where Conv denotes the closed convex hull. If Hq is normal and A is normal for 
every A & M. then 

Spec{H) C Conv(E) + Conv(Spec(M)) (1) 

Proof Theorem 1 implies that for each A E M. 

E + Spec(^) = Spec(//o ® / + / ® ^) C Spec(i/) 

and this yields the first inclusion. The second depends on use of the numerical 
range to give 

Spec(i/) C Num(i7) 

C Num(/7o) + Num(l^). 

Now z lies in the numerical range of V if and only if there exists / G 1^{X,}C) 
of norm 1 such that z = {Vf,f). Putting = /(a;)/||/(a;)||, provided this is 
non-zero, and iJ,x = \\f{x)\\^, we see that ^ is a probability measure on X and that 

z= l^x{Vxgx,gx) e Conv(Num(M)). 

xex 
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Hence Num(y) C Conv(Num(M)), and the first statement of the theorem fol- 
lows. The second statement is a consequence of the fact that Num(S) equals 

Conv(Spcc(i?)) for any normal operator B. 

Let B{x, r) denote the closed ball {y : — 1/| < r}. The next theorem complements 
Theorem 3. 

Theorem 4 If A is normal for every A & M. then the spectrum of H satisfies 

Spec(/J) C Spec(M) + 5(0, e) (2) 
where e = ||-ffo||- If Hq is normal then 

Spec(//) <ZE + B{0, li) (3) 
where n — max{\\A\\ : A G Ai}. 

Proof If V is normal then using Spec(\^) C Spec(A4) we see that 

11(^-^7)-^ =dist(^,Spec(y))-^ <dist(z,Spec(M))-^ 

for all z ^ M. Since z ^ M + -8(0, e) is equivalent to dist{2;, Spec(A4)} > ||i7o||, 
it implies 

||//o||||(^-^/)-^||<l 

and the resolvent expansion for (Hq + V — zl)~^ is norm convergent. The proof of 
the second part of the theorem is similar. 

We also wish to classify the spectrum of nsa operators acting on P{X, JC), and for 
this purpose we assume that X is provided with a metric d such that every ball 
B{x,r) = {y & X : d{x,y) < r)} is finite and such that T acts as a group of 
isometrics of X. Given a function f : X ^ K, with ||/||2 = 1 we define its variance 

by 

var(/) = nii^ ^ d{x,yf\f{x)\^ 

^ x&X 

and its expectation to be any of the points in X at which the minimum is achieved. 
The following theorems have analogues in which the variance is replaced by higher 
order moments, or suitable subexponential weights. 

Lemma 5 // 1|/||2 = 1 and 

v{x)^Y.d{x.yf\f{y)? 

y&X 

is finite for some x E X then it is finite for every x E X and v{x) increases 
indefinitely as x ^ oo. Thus the minimum of v{-) is achieved at a finite number 
of points only. If Xi, i = 1,2, are points at which v has the same minimum value 
s then d{xi,X2) < 2s^/^. 



5 



Proof If f (x) < oo then for any -u G X we have 



v{u) < 2 Y.{di^^yf + d{x,uf]\f{y)\^ = 2{v{x) + d{x,uf} < oo. 

y&X 

by the triangle inequahty. If the finite set F satisfies 

y&F ^ 

then 

v[x)>Y.d[x,yf\f[y)\^>\d{x,Ff 

y&F ^ 

which increases indefinitely as x — > oo because of our assumption that all balls of 
finite radius contain only a finite number of points. 

Now suppose that s = min{t>(x) : x G X} and that v{xi) = f (X2) = s. Then by 
the triangle inequality 

2s = Y.{di^^^yf + dix2,yr}\f{y)\' 

> Wd{xuX2r\fiy)\' 

= -d{xi,X2y 
which implies the second statement of the lemma. 

Following p we define the localized spectrum aioc{A) of any bounded operator A on 
/^(X, /C) to be the set of all A G C such that there exists a sequence /„ G /^(X, /C) of 
unit vectors such that \\Afn — A/„|| — >• while var(/„) remains uniformly bounded. 
If A is an eigenvalue then one would expect its corresponding eigenfunction to 
decrease rapidly at infinity and hence to have finite variance, in which case A 
would lie in a\oc{A). What is more surprising is that o"ioc(^) can be much larger 
than the set of eigenvalues of A. 

Theorem 6 If H = Hq + V where V is (F, Ai) pseudo-ergodic and K = Hq + W 
where W : X ^ M. is arbitrary, then 

C"loc 

Thus every eigenvalue of K lies in the localized spectrum of H . Moreover ifV,W 
are both {r,Ai) pseudo-ergodic then they have the same localized spectrum. 

Proof First note that if / G P{X,}C) has unit norm and 7 G F then g = f/y/ 
has the same variance as / because F acts as a group of isometries of X. It is 
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a consequence of the definition of pseudo-ergodicity that there exists a sequence 
7(n) G r such that Hn = U~^l^-^HU^(^n) converges strongly to K. Now let ||/m|| = 1, 
var(/m) < s and \\Kfm — A/m|| < for all m G Z+. Given m 

||-f^(^7(n)/m) - A(t/^(n)/m)|| 



as n oo. Therefore there exists n(m) such that = f^7{n(m))/m satisfies 

m 

for all m G Z+. Since varlgm) < s for all m it follows that A G aioc{H). 

We next turn to the essential spectrum. We say that z lies in the essential spectrum 
of a bounded operator A if A — zl is not a Fredholm operator. We will need the 
following known result. 

Proposition 7 Suppose that z G C and for all e > and all finite N there exists 
an orthonormal set fi, such that \\Afn — zfn\\ < s for all 1 < n < N . Then 
z lies in the essential spectrum of A. 

Proof Suppose that z G C satisfies the conditions of the proposition. If ker(y4 — 2;/) 
is infinite dimensional then A — zl is obviously not Fredholm, so let dim(ker(74 — 
zl)) < N where is finite. The assumption implies that for all e > there exists 
an A^-dimensional subspace L such that f E L implies 

P/-^/ll< ^11/11- (4) 

Because dim(L) > dim(ker(A — zl)) there exists / _L ker(74 — zl) such that 
holds. Since e > is arbitrary, A — zl cannot be Fredholm. 

Lemma 8 Suppose that there exists a (r,A^) pseudo-ergodic potential V on X 
where Ai C £(/C) contains more than one point. Then for any finite subset F 
of X and any finite N there exist 71, ...,7Ar G F such that {7n-F}^=i are pairwise 
disjoint. 

Proof Let us first put N = 2. Let mi,m2 G Ai and \\mi — m2\\ = 25 > 0. Also let 
W : F ^ C{IC) satisfy W{x) = rrii for all x E F. Since V is (F, Ml) pseudo-ergodic 
there exist 7^ G F such that ||l^(7ja;) — < 6 for all x G F, or equivalently 
II 1^(1/) — mill < S for all y G •jiF. This implies that 71F fl 72 -F = 0. 

We next prove that if the lemma holds for then it holds for 2N; we can then 
complete the proof by the use of induction. We put F = [jjLi IjF and let /?2 ^ P 



l^7(n)-^^7(")/"i A/„ 
\Hin fm A^mll 



11^/. 



A/m|| < — 

m 
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be such that PiF fl = 0. This yields the statement of the lemma for the sets 
Pi'jjF where i = 1,2 and 1 < j < A^- 

Theorem 9 If H = Hq + V where V is (T,M.) pseudo-ergodic and M. contains 
more than one point, then H has no inessential spectrum. 

Proof If A G Spec(if) then either (i) for every £ > there exists / G /^(X, /C) such 
that 11/11 = 1 and \\Hf - A/|| < e, or (ii) for every e > there exists / G /^(X, /C) 
such that 11/11 = 1 and \\H* f — A/|| < e. We assume (i), the proof for (ii) being 
similar. By approximation we may assume that each / has finite support F. Now 
for any £ > and any finite N let 71, ...,7Ar G F be such that 7jF are pairwise 
disjoint. Put F = U^Ii 7rF and define 1^ : F ^ M by W{-fix) = V{x) for all 
X E F. Since V is (F, A^) pseudo-ergodic there exists 7 G F such that 

\\Vijy)-Wiy)\\<e 

for all y E F. Thus 

\\V{^^,x)-V{x)\\<e (5) 

for all X G F and I <i <N. 

We now put fi{x) = f{%^l~^x) for all x G X and observe that /j have supports 
within '~f'~fiF, which are disjoint, so {fi}^i form an orthonormal set. It follows 
from condition (i) and (^) that 

\\Hf,-XM\<2e 

for all 1 < i < X. This implies that A lies in the essential spectrum of H by 
Proposition 7. 



3 The nsa Anderson model 

In this section we apply the above ideas to an example of physical and biological 
importance. We first consider the one-dimensional nsa Anderson operator 

Hfn = e-^/„-l + eVn+l + Vnfn (6) 

acting on /^(Z) (so that /C = C), where g > and V is a. (F, M) pseudo-ergodic 
potential, F being the group of all translations of Z and M being a compact subset 
of C. The potential V may be generated by assuming that its values at different 
points are independent and identically distributed according to a probability law 
which has compact support M. 

Fourier analysis quickly establishes that Hq is normal with spectrum the ellipse 

E = {e3+'^ + e-'-'^ : ^ G [0, 27r]} (7) 
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following which Theorem 3 implies that 

E + M C Spec(if) C Conv(E) + Conv(M). (8) 

A more precise determination of Spec (if) depends upon the size of g, the choice 
of M and the use of Theorem 6, extending what we already proved in [Q. Given 
any finite sequence a = (ao, ai, . . . , an-i) £ M'^ let Wa be the periodic potential 
such that Wa,m = ar ii m = r mod n. The eigenvalue equation 

e-Vm-l + Wc,,mfm + fm+1 = Xfm (9) 

may be rewritten in terms of Wm = (/m-i, fm) ^ as Wm+i = WmAm where 







1 e-9{X-Wc, 



Thus 



Wn(r+1) = WnrB 

for all r G Z where B is the transfer matrix 



B = AqAi . . . An-l. 



Since 



det(5) = n det(A.) = e^ 



-2ng 



r=0 



it follows that at least one of the two eigenvalues /ii,/i2 of B satisfies < 1. If 
we write 

bu{X) 612(A) 
&2i(A) 622(A) 



B 



then one may prove by induction that &22(A) is a polynomial of degree in A while 
the other coefficients are of lower degree. 

The solution / of (^ corresponding to an eigenvalue fi of B is exponentially in- 
creasing or decreasing on Z according to whether > 1 or < 1 respectively. 

Theorem 10 Let E"^ denote the ellipse 

= {e*^ + e-^'^a-ie . q ^ [_^^ 

and let 

E„ = {A :6n(A)+ 622(A) GE"}. (10) 

Then B has an eigenvalue of modulus 1 if an only if X G Ea- Moreover E^ is 
closed and hounded with 

Ea C Spec{H). 
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Proof If /ii = e*^ for some 9 G [— vr, vr] then /X2 = e and 

6n(A) + 622(A) =e^^ + e-2"3-^^ 

or equivalently A G -Eq. The converse also holds. Our comments above on the 
degrees of hij{\) imply that + /i2| increases indefinitely as |A| grows. Therefore 
one of the /Xj must have modulus greater than 1 for large enough |A| and such A 
cannot lie in Ea-, therefore must be bounded. The fact that E^ is closed follows 
directly from its definition. 

Corresponding to any \ E Ea there exists a solution / of (H) such that fm+n = fm 
for some 6' G R and all m G Z. This / is bounded but its norm is infinite. If we 
put 

f _ p-£|m| f 
J e,m ^ Jm 

then a direct and well-known calculation shows that ||/e||2 00 and 

11(^0 + iy)/,-A/,||2 
II/.II2 

as £ ^ 0. Applying Theorem 1 we deduce that 

A G Spec(ifo + W)C Spec{H). 

The set C \Ea is the union of disjoint components and the number of eigenvalues 
fij of B which have modulus less than 1 cannot change within each component, 
because the eigenvalues depend continuously on A. This number must be either 1 
or 2, and within the unbounded component it is 1. The following theorem joins 
the components into two sets. 

Theorem 11 If X lies in 

/a = {A:6n(A) + 622(A)Gint(i?")} (11) 

then all solutions of ^ are exponentially decreasing. If, however, A lies in 

Oa = {\: 611(A) + 622(A) G ext(E")} (12) 

then there is an exponentially increasing solution of ([^. The three sets la, Oa and 
Ea are disjoint and cover C. 



Proof The condition ([TT|) holds if and only if both /ij have modulus less than 1, and 
this implies that every solution of (^) is exponentially decreasing on Z. Similarly 
The condition (|12[) holds if and only if one fii has modulus greater than 1, and this 



implies that one non-zero solution of (P) is exponentially increasing on Z. 
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The explicit description of the above sets depends upon the value of n. For n = 1 
we have a G M and 

Ea = E + a. 
If n = 2 and a = (cto, ai) G then 

^ ^ [ e-a{X - ao) e'^^^A - ao)(A - ai) - 1} 

and Ea is the set of A such that 

e-2^{(A-ao)(A-ai) -2} G E'. (13) 

This equation may be solved to present A explicitly as a function of Q. For larger 
values of n it is probably only practicable to find E^ numerically. 

The special case = p = 1 of the following theorem was proved in The idea 
owes much to the theory of block Toeplitz matrices [0, 0, |, 0. 

Theorem 12 Let H he defined by ^ where g > and V is a (Z, M) pseudo- 
ergodic potential. If a E M" and (3 G then 

Proof We consider the operator K = Hq + W acting on /^(Z) where 




ar if m > and m = r mod n 
(3r if m < and m = r mod p. 



We then consider the solutions of 

where A G fl Op. Since A G Op there exists a solution / which is exponentially 
growing for m < 0, i.e. which decreases exponentially as m — >■ —oo. Continuing 
this solution to positive m it follows from A G that / also decreases exponentially 
as m ^ oo. Hence / is an eigenvector of finite variance and A G crioc(-^) ^ <y\oc{H)- 

Theorem 13 // in addition to the hypotheses of the last theorem we put M = 
[—fi,fi] then Spec(if) = E + [— /i, /i] for all /i > e^ + e~^. Moreover G Spec(if) 
if and only if fi > — . 

Proof The first statement only needs the observation that the two sides of ® 
coincide under the given condition. If fi < — then ^ Spec (if) by Theorem 
4. Now G if and only if e~^'^(— /x^ — 2) G E"^ by (p^ , and this is equivalent 

to /X = e^ — e~^; for such fj, one has G Spec (if) by Theorem 10. For smaller fi we 
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have G and for larger /i we have G 0(_^,^). Therefore G /(o) H 0(_^,^) 

for yU > e^ — e~^, and G crioc{H) by Theorem 12. 

If M = [— /i, /i] the above theorems admit the possibihty that there are two holes 
in the spectrum on either side of the origin for 

e^ - e"^ < /i < e^ + e"^. 

We nevertheless conjecture that one has Spec(iJ) = Conv(£' + M) for all /i > 
e^ - e-5. 

We contrast the above with the case in which M = {±/i}. The following theorem 
completely determines the real part of Spec (if) under the stated conditions. 

Theorem 14 // M = {±fi} and /i > e^ + e~^ then 

(Conv(E) + iJ,)U (Conv(E) - jJ,) C Spec{H) C B{fi, e^ + e~») U e^ + e'^) 

and 

Spec(i7) C Conv(^) + 

Proof The first inclusion of the statement follows from the case n = p = 1 of 
Theorem 12 as in [Q. The second follows from the first half of Theorem 4, and the 
final one follows from Theorem 3. 

We conjecture that the first inclusion is actually an equality. 
Corollary 15 If M = {±yu} then G Spec(ii') if and only if 

- e'^ < 12 < + e"^. 

Proof If /i < e^ — e~^ then ^ Spec(ii) by combining Corollary 2 and Theorem 
13. If /i > e^ + e-3 then ^ Spec(if) by Theorem 14. If /i = e^ - e'^ then 
G ^(-;.,^) C Spec(ii) by Theorem 10. If /x = e^ + e'^ then e C Spec{H) by 
Theorem' 10. Finally if - e'^ < fi< - e'^ then G 0(_^,^) n C Spec(ii) 
by Theorem 12. 

We next turn to the nsa Anderson model in Z". The operator H on /^(Z") is 
defined by 

{Hf){m, n) = {Hof){m, n) + V(m, n)f{m, n) 

where 

{HQf){m, n) = f{m + 1, + f{m — l,n) + f{m, n + 1) + f{m, n — 1) 

for some g > 0. We assume that V is real- valued and pseudo-ergodic with values 
in M = [— It follows by Fourier transform methods that Hq is normal with 
spectrum equal to 

E = E+[-2(n-l),2(n-l)] 
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where E is the elhpse defined by (|^). This set is connected with a hole around the 
origin if n = 2 but it may or may not have such a hole for n > 3. This phenomenon 
is a result of the particular choice of lattice used to discretize the Laplacian. If fi 
is sufficiently small the same applies to Spec (if). 

Theorem 16 If fi > + e^^ — 2{n — 1) then Spec(if) is the convex set 

E + [-(x- 2(n - 1), /i + 2(n - 1)]. 

Proof As in the one-dimensional case we need only observe that the two sides of 
(H) are equal under the hypotheses. 

We next mention the same operator acting in /^(X) where 

X = {{m,n) : m e Z,l < n < N} 

subject to Dirichlet boundary conditions; the Neumann case is similar. We may 
carry out an analysis similar to that above if we are only concerned to determine the 
spectrum, but more detailed spectral information is obtained by putting /^(X) = 
/2(Z, /C) where /C = C^. We then put 

{Hof){m) = e'f{m + 1) + e'Vl^ - 1) 

and 

V{m){r, s) = 




where 1 < r, s < in all cases. Note that Hq is normal and V{m) is a self-adjoint 
matrix for all m G Z, so all of the theorems of Section 2 apply. Using such ideas it 
is possible to analyze the localized spectrum of H as in the one- dimensional case. 

We finally comment that certain random bidiagonal operators can also be treated 
by the methods of this paper by making the appropriate choice of Hq, as can a 
variety of other operators whose matrix coefficients depend only on m — n whenever 
m ^ n. See [||, & |l8i, which use probabilistic rather than pseudo-ergodic methods. 



4 Resolvent Norms 

The spectral behaviour of a bounded operator A acting on a Hilbert space Ti can 
be measured in several ways. In pseudospectral theory one examines the contours 
of the function 

\\{A-z)-^\\-^ if2^Spec(A) 
ifzeSpec(A). 
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This function converges to zero as z approaches the spectrum of A because of the 
upper bound 

s(y4, z) < dist(z, Spec(A)) 

and the case of most interest is when s{A, z) is very small for z far from the 
spectrum. The determination of the pseudospectra, defined as the family of sets 
{z : s{z) < e} for all positive e, is computationally heavy, but the family carries 
much more information than the spectrum alone |jl], ^, |l^, |15|, |16|, p!7| . 

Lemma 17 The function s{A, ■) satisfies the Lipschitz inequality 

\s{A, z) — w)\ < \z — w\ 

for all z,w E C. 

The proof uses the formula 

.(A^)=mf{||(A-^)/||/||/||:0^/G7^} (14) 

valid for all z ^ Spec (A). Note that this may be false for z G Spec (A), as one may 
see by considering the operator A on /^(Z"*") defined by 



Af{n) 



if n = 1 

f{n-l) ifn>2. 



The next theorem provides an upper bound on s{A, ■ ) which may be used to 
compute it numerically. Let L be a finite-dimensional subspace of Ti and let P be 
the orthogonal projection onto L. We define B{A, L, z) to be the restriction of 

P{A - ziyP{A - zI)P + PA*{I - P)AP 

to the subspace L, and a{A, L, z) to be the square root of the smallest eigenvalue 
o{B{A,L,z). 

Theorem 18 // we put 

s{A, L, z) = min{cr(74, L, z), <j{A*, L, z)} 

then 

s{A,L,z) > s{A,z). 

The functions s{A, L, ■ ) decrease monotonically and locally uniformly to s{A, ■ ) 
as the subspaces increase. 



14 



Proof It follows from its definition that 



<7(A L, z) = min{||(A - ;^)/||/||/|| : 7^ / G L}. 
It is clear from this that (t{A, L, z) decreases monotonically and pointwise to 

a{A,z)=mi{\\iA-z)f\\/\\f\\:Oy^fen}. 

If z ^ Spec{A) this equals s{A,z). Similar comments apply with A replaced by 
A*, and we also have (t{A, z) = a{A*,z) for all z. 

On the other hand if z E Spec(A) we have either a{A,z) = or (t{A*,z) = 0, 
or both. This implies that s{A, L, z) converges monotonically and pointwise to 0. 
Since all the functions involved are Lipschitz continuous with Lipschitz constant 
1, the convergence must be locally uniform. 

Now suppose that H equals P{X, K) and L is defined as the space of all functions 
with support in a particular finite region fi. The above theorem is better than 

the mere computation of the spectrum of PAP restricted to L (possibly subject 
to certain boundary conditions on dO) because it gives rigorous upper bounds to 
s{A, z) rather than uncontrolled approximations. Another advantage is that it 
provides an upper bound for 5(^4, z) for every extension of the operator A beyond 
the subspace L. Because of its approximate nature one cannot determine the 
spectrum of A exactly using the above theorem, but it may be possible to get good 
approximations to the pseudospectra, which are often of greater importance for 
such operators. 

We now turn to pseudo-ergodic operators, working in the technical context of 
Section 2. The following theorem indicates how one may get rigorous upper bounds 
and approximations to the pseudospectra by selecting appropriate potentials W . 

Theorem 19 If H = Hq + V where V is (F, M) pseudo-ergodic and K = Hq + W 
where W : X ^ M is arbitrary, then 

\\{H-zI)-'\\>\\{K-zI)-'\\ 

for all z & C. Therefore 

s{H, z) = mm{s{Ho + W,z):W E M^}. 

IfV,W are both (F, M) pseudo-ergodic then the resolvent norms and hence pseu- 
dospectra of H and K are equal. 

Proof By Theorem 1 we need only consider the case in which z does not lie in the 
spectrum of either operator. If s{K, z) < c then there exists / e t^iX, K) such that 
IKX — z)/|| < c||/|| and by approximation we may assume that / has finite support. 
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Using the pseudo-ergodic property of H there exists g of finite support such that 
\\{H — z)g\\ < c\\g\\ and this imphes that s{H,z) < c. Hence s{H,z) < s{K,z). 
The remainder of the proof follows Theorem 1 or Corollary 2. 

For the nsa periodic Anderson model with M = [— /i, fi] the asymptotic limit of 



the finite volume spectrum has been determined |T0[, and it is seen that for certain 
ranges of the parameter fi zero does not lie in the asymptotic spectrum, which is 
the union of a set of complex curves. On the other hand the spectrum of the same 
operator on any finite interval subject to Dirichlet boundary conditions is entirely 
real. It has been suggested in that for periodic boundary conditions there 
is no pseudospectral pathology of the type which occurs for Dirichlet boundary 
conditions. However, our results demonstrate that spatially rare special sections of 
a random potential have a dominant effect on the spectrum of the infinite volume 
nsa Anderson operator. This should not be taken as an indication that our results 
are unphysical: it is well known that the behaviour of bulk materials is often 
radically affected by the presence of low concentrations of impurities and/or defects, 
and one should expect the mathematics to reflect this. 

We have implemented the above ideas numerically using Matlab for the operator 
H defined by where e^ = 2 and Vn are independent random variables uniformly 
distributed on [—3,3]. We took L to be the subspace of all sequences with sup- 
port in [1, 100] and computed the minimum value of ct{H, L, x) over 1000 different 
choices of the potential V. We chose to study real x G [0,6], but complex values 
of X in any region can be accommodated by the same method. This yielded the 
upper bounds s{H, x) as follows (the omitted values of s{H, x) all vanish to the 
given accuracy). 



X 


s{H, x) 


0.0 


0.0283 


0.5 


0.0203 


1.0 


0.0084 


1.5 


0.0015 


4.5 


0.0044 


5.0 


0.2233 


5.5 


0.6259 


6.0 


1.0817 



Our general theory shows that the real part of the spectrum of this operator is 
[—5.5, 5.5], which is consistent with the numerical conclusion that 

\\{H-xI)-^\\ > 10^ 

for all 1.0 < X < 4.5 and 

\\{H-xI)-^\\ > 10^ 
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for all 2.0 < x < 4.0. (Of course the numerical calculation can also be carried out 
in cases in which one does not have a prior theoretical solution!) The eigenvectors 
of B{A,L,x) corresponding to the smallest eigenvalues were also computed for 
several values of x. As expected from the theory of localized spectrum, they were 
all highly concentrated around some point in the interior of [1, 100], and negligible 
at the ends of the interval. 

We finally examine the behaviour of the resolvent norm at the point z = 0. To be 
precise we consider the nsa Anderson model with M = [— /i, /i] acting on /^(Z) for 
various values of fi. Recall that Theorem 13 states that G Spec{H) if and only if 
/i > e» - e-3. 

Theorem 20 // < ii< - e'^ then 

\\H-^\\-^ = - e-^ - X. 

Proof If we exhibit the fi dependence of H explicitly and put = s{H^, 0) then 
it follows from (p^ that 

|t(yU) — t(z/)| < l/U — 

for any < yU, z/ < e^ — e~^. Since t{0) = e^ — e~^ and t{e^ — e"^) = we conclude 
that 

t{fx) =q9 _Q-9 _ ^ 

for all < /i < e^ -e-5. 

5 Constrained Potentials 

We have avoided the use of any probabilistic methods by the introduction of the 
concept of pseudo-ergodicity. We now explore the variety of situations in which our 
ideas are applicable. The obvious possibility is to assume that is a probability 
measure with support equal to the set M. C £(/C) and to assume that are inde- 
pendent random variables as a; G X varies and that each is distributed according 
to yU. However, even if we assume that are independent, we may permit each 
Vx to be distributed according to a different probability measure with support 
equal to A^. These measures need not even be F-stationary, but they must satisfy 
the following condition. For every open set U C C{JC) such that U f\ M. ^ ^ there 
must exist a constant cu > Q such that [ixiU] > cu for all x G X. This is sufficient 
to imply that V is (F, A^) pseudo-ergodic almost surely by the usual probabilistic 
argument. For all such probabilistic models the spectrum (or localized spectrum) 
of the operator H is the same. 

Similar remarks apply to a variety of other probabilistic models in which the values 
Vx are not independent. There is one situation, however, in which changes in 
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the spectrum may arise. We say that a potential V satisfies the local constraints 
Q = (A^, 7i, 7fc, A/i, A/fc) where 7i G F and A4,Mi are closed, bounded subsets 
of C{}C) under the following conditions. For all a; G X we require that Vx E 
and also that 



for all i = 1, k. Even more general constraints can be formulated. We then say 
that V is (F, Q) pseudo-ergodic if it satisfies the constraints Q and for any other 
potential W which satisfies the same constraints and any finite subset F of X and 
any e > there exists 7 G F such that 

11^7^ - Wx\\ < e 

for all X E F. These constraints force a relationship between the values of Vx at 
neighbouring points which is stronger than a mere probabilistic correlation. 

Lemma 21 If Hj = Hq + Vj where Vi is (F, A^) pseudo-ergodic and is (F, Q) 
pseudo-ergodic, then 



Any two (F, Q) pseudo-ergodic operators have the same spectrum. 

Proof The first statement is a consequence of Theorem 1. The second involves 
adapting the proof of the same theorem. 

We now apply the above ideas in a simple context. We assume that X = Z, that 
F is the usual translation group acting on Z, and that /C = C. We assume that 
a, h are two positive constants and impose attractive constraints Qi of the form 



for all G Z. Although we are not able to prove Theorem 12 in full generality 
under such conditions the important special case n = p = 1 is still valid. 

Theorem 22 Let H he defined by ^ where g > and V is a (Z, Qi) pseudo- 
ergodic potential. We have 

E + [-a, a] C Spec(iJ) C Conv(E) + [-a, a] 

where E is given by Ifa,/3 G [—a, a] then 

Proof The first statement of the theorem is proved as in Theorem 3. For the second 
part we follow the method of Theorem 12 but for the operator K = Hq + W acting 
on /^(Z) where 



Vx - V^^x G Mi 



—a < Vn < a 



Vn - Vn+l\ < b 
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Here we take large enough to ensure that W satisfies the constraints Qi. 

A more interesting variation upon our earher theory occurs if we impose the re- 
pulsive constraint Q2 defined by 

-a<Vn<a , \Vn-Vn+i\>b 

for all n G Z, where < 6 < 2a. This excludes constant potentials, thus rendering 
the first inclusion of Theorem 3 invalid. The range of a (F, Q2) pseudo-ergodic 
potential V is equal to M = [—a, a — b]U [b — a,a]. 

The spectrum of the Anderson model (|) is easy to determine in the self-adjoint 
case, and we start with this. 

Theorem 23 If g = and V is (F, Q2) pseudo-ergodic then the spectrum of the 
operator H = Hq + V defined by (Qj is given by 

Spec{H) = T U (-T) 

where 

T = 



b 62 



Thus Spec(if) has a spectral gap if and only if a < b < 2a. 

Proof Let W be the potential Wn = (— 1)"&/2, so that W + si satisfies the con- 
straints Q2 for all real s such that |s| < a — 6/2. It follows from Theorem 1 that 
ifS = Spec{Ho + W) then 

S +[b/2-a,a-b/2]CSpec{H). (15) 

Conversely \\V — W\\ < a — b/2, so the perturbation theoretic argument used in 
Theorem 4 implies that 

Spec(if) CS+ [b/2 -a,a- b/2]. 

We deduce that 

Spec{H) = S+ [b/2 -a, a- b/2] 

and complete the proof by using a Bloch wave analysis to compute the set S. 

Now let us denote the same operator by Lg for g > 0. We may regard Lg as a 
perturbation of Lq and use the argument of Theorem 4 to show that 

Spec(Lg) C Spec(Lo) + 5(0, e^ - 1). 

We may also use Theorem 4 as it stands to obtain an outer estimate of Spec(Lg). 
We may obtain inner estimates by the method of Section 3 provided we are careful 
to avoid the use of constant potentials. 
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Theorem 24 We have 




— a,a C Spec(Lg) 



where 



S 




+ 2 + e29+»« + e-2s-»e : 9 e [-n, n] 



Proof if we put oq 



6/2 and ai = 6/2 and solve (0) for A we obtain 

E(-b/2,b/2) = S. 



The remainder of the proof follows Theorem 23, using the last part of Theorem 10. 

Note that for small positive g, S consists of two closed curves on opposite sides of 
the y-axis, but for large 5^ it is a single curve enclosing the origin. 
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